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OO ; Abstract 

Q ■ We develop exact field theoretic methods to treat turbulence when the effect of 

pressure is negligible. We find explicit forms of certain probability distributions, 
demonstrate that the breakdown of Galilean invariance is responsible for intermit- 
tency and establish the operator product expansion. We also indicate how the effects 
Q^' of pressure can be turned on perturbatively. 
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Turbulence is an old and tantalizing subject. Enormous amounts of data and ideas 
have been accumulated during this century and still the problem is not solved. In our 
opinion, the reason lies in the fact that the necessary field-theoretic tools have appeared 
only recently. 

Two years ago an attempt was made to apply the methods of conformal field theory 
to the case of two-dimensional turbulence [1]. The main concepts of this work were 
the following. First, one looks at the steady state condition, which relates the equal 
time, N- point and N-l-1 - point functions. Then one argues that in the inertial range 
these relations can be solved exactly by field theories satisfying fusion rules or operator 
product expansions (OPE). There appeared to be infinitely many solutions. An additional 
constraint on these solutions follows from the constant flux conditions. 

It has been noticed (although not really exploited) in [1] that there exists a striking 
analogy between the constant flux states in turbulence and axial (and other) anomalies 
in quantum field theory. The latter are violations of the naive conservation laws caused 
by the ultra-violet regularization. In the case of turbulence the ultraviolet regularization 
arises from viscosity and results in an energy flux through the inertial range. When 
the steady state condition with these two tools was analysed it appeared that the third 
ingredient was needed. Namely the physical correlation functions contained so-called 
condensate terms, which were S - functions in the momentum space and represented the 
large scale motions of the fluid. Their role was to cancel infrared divergencies yhat arose 
from the fleld-theoretic fluctuations. Precise form of these terms depends on the large 
scale region were the energy is pumped into the system. The task of joining the inertial 
range with this region remained unsolved in [1] . 

It is highly desirable to have an exactly soluble model in which the above ideas can be 
tested at work. In this article we will discuss such a model, which also is of independent 
physical interest. The model in question is simply the Navier-Stokes equation with white 
noise random force and with the pressure set equal to zero. In one dimension this is 
known as the Burgers equations. 

Such equations have been exploited in the past in many different physical situations 
(like galaxy formations [ 2 ], crystal growth [ 3 ] etc). Recently they were the subject of 
deep mathematical investigations [ 4 ] 

In a remarkable paper [ 5] some striking numerical data concerning Burgers turbulence 
were obtained and an appealing qualitative picture of the phenomenon has been proposed. 
This work to a large extent inspired my interest in turbulence without pressure. Another 
important work in this area is the recent paper [6] on which I shall comment later. 

In the present paper we shall formulate a general new method for analyzing the in- 
ertial range correlation functions, based on the ingredients mentioned above (OPE and 
anomalies). The method with minor modifications is also applicable to the problem of 
advection of passive scalars and other cases. It is obvious that the ideas we develop below 
will become a part of the general theory of turbulence. They may also have a considerable 
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"back reaction" on the field theory. 

Let us start with the one dimensional case. Burgers equation has the form: 

Ut + UU^ = Z/M^^ + f{xt) (1) 

< f{x, t)f{x, t) = k{x — x)5{t — t) 

Here the function k defines the spatial correlation of the random forces. Consider the 
following generating functional: 

Z{\iXi \ . . . \\nXn) = (exp^ Ajti(a;j^)^ (2) 

From (1) we derive : 

Z + Y.^3Q^iY,d^) = Y. \ (ifi^ji) + ^«") exp E >^ku{x,t)) (3) 

The first term in the right hand side is easy to treat since the force f{xt) is Gaussian and 
we can use the standard trick [ 7] of the theory of Langevin equations: 

(^f{xt) exp E -^j^l^iO) — X/ '^(•^ ~ (4) 
Our equation takes the form: 

F) 1 87 
By D we denote the dissipation term: 



D = z^E \ (u"{xjt)expY^ Xku{xkt)^ (6) 

If the viscosity were zero our task would be completed since we have a closed differential 
equation for Z. To reach the inertial range we must, however, keep u infinitesimal but 
non-zero. The anomaly mechanism mentioned above implies that infinitesimal viscosity 
produces a finite effect, whose computation is one of our main objectives. First, however, 
let us transform and interpret the inviscid equations (5) ( dropping the D-term). 

Let us introduce the function F given by : 



Z — Ai . . . X]\fF(^XiXi . . . Ajv^^Af) 



(7) 



We have: 



^ + E^^-E-(-.-.)W = ^ 



(8) 



Here: 

b ^ 



We can now introduce the Fourier- transform JF = J^{uiXi . . . uiyx^) 
which satisfies: 

obtained from(8) by the substitution A ^ ^ 

The function T has a simple meaning. It can be interpreted as: 

J^^{0{uiXi)...e{uNXN)) (10) 

where 9{ukXk) — 9{uk — u{xkt)) 

and the last 9 being a step function. In order to get the N-point probability distribution 
Z one has to differentiate according to (7): 



QN 

Zn = j; ^^iv (11) 

OUi . . . OUn 

We supplied here our correlation functions with the subscript A'^ to indicate the number 
of points on which these functions depend. 

Of course, equation (9) could have been obtained directly by computing time derivative 
of the 9 - field. It is also easy to write to express the Dj^ - term through J-n+i By using 
(6) we obtain: 

^ ^'^J uoduo-^J^N+i{uo,Xj + y;ui,xi . . .un,xn) \y^o (12) 

Equations (9) and (12) give a chain of relations remarkably similar to the BBGKY 
equations of statistical mechanics [7]. One can hardly hope to solve these equations 
exactly. But we are interested in the inertial range, which means that we have to take 
the limit — > 0. We will show now that in this case the system of equations closes and 
gives us an equation for turbulent kinetics, much in a same way as the Boltzman equation 
becomes exact in the limit of small densities. 

The main ideas of the derivation are the following. For the large Reynolds numbers, 
corresponding to small v there are two relevant scales. The first, L , is defined by the 



3 



size of the system and provides an infrared cut-off. The second ,a ^ L, is the scale 
at which dissipation becomes relevant. The ratio ^ goes to infinity together with the 
Reynolds number. By the existence of the inertial range we mean the conjecture that 
the correlation functions have a finite limit at zero viscosity, provided that we keep 
Xi — Xj fixed. They can have singularities at coinciding points, which must be understood 
as being smeared by the viscosity at the scale a . In fact, this scale is determined by 
the condition that as we let v go to zero, dissipation remains finite. This means, that 
we have to find the leading singularity in (12) as y — > 0, and compensate for it by an 
appropriate scaling of z/(|/) . All sublcading terms will give vanishing contribution in the 
inertial range (in the limit of the infinite Reynolds numbers). The task of finding the 
leading singularities is precisely what the OPE was developed for. 

However, we should warn the reader that what follows is essentially a self - consistent 
conjecture. In the case of statistical mechanics, when deriving hydrodynamics from the 
BBGKY equations it is necessary to assume the decrease of the correlations [7 ], a self- 
consistent assumption that is difficult to prove from first principles. In our case this 
property is replaced by the OPE. 

To understand how they work, let us reexamine the derivation of the previous equa- 
tions. They were based on the fact that modulo the stirring force and the viscosity we 
have a sequence of conservation laws: 



^(M") + ^^(«"+^)RiO (13) 

(The sign ^ here means that we don't write terms coming from the viscosity and the 
stirring force) 

Equations (5) and (9) can be interpreted as relations for the generating functionals 
{u"'^{xi) . . (xfc)). They involve both the stirring force and the viscosity. The former 
was already accounted for, while the latter presents a problem. The main rule of the game 
is that in any equation , involving space points separated by the distance larger than a 
viscosity can be set to zero. Thus it is perfectly legitimate to use the inviscid limit for 
n=l, because in this case we exploit the steady state condition: 

^{u{xi) . ..u{xn)) = 

\ Xi — Xj \':$> a (14) 

However, starting from n=2 we have a problem, since in this case we have to take time 
derivatives of the product of w - s at the same point. To circumvent this problem in the 
case of n=2 let us replace: 

u'^(x) =^ m(x + |)ii(x - |) 

\ Xi - Xj |> y > a (15) 
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and let y — after the viscosity is taken to zero. In this case the use of the inviscid 
equations is justified, but we will get an anomaly in the conservation law, due to the 
point splitting. We have after simple algebra: 

- j^H^ + - |)) - ^^(^'(^1)^(^2)) + (1 ^ 2) « + ""o^"^)' (^^^ 

Xi^2 = X±- 

Here we have introduced the first dissipative anomaly operator: 

1 d 1 1 

aoix) = lim - — (u(x -\ — y) — u(x y))^ (17) 

In deriving this formula we set y to zero inside all terms containing x- derivative. This is 
possible because all the correlation functions have finite x -dependent limit at zero y. We 
also used the identity: 

|(«'(. + f ) - - |)) = + |) + _ |„ ^ 

The anomaly would be zero if u{x) were differentiable, since then the RHS of (17) is 
~ y"^. However, the steady state condition dictates the opposite. Indeed, one of the 
consequences of eq. (5) is that in the steady state: 

|(«2) = «:(0)-(ao) = (18) 
and hence we have the famous Kolmogorov relation: 

[u{xi) - u{x2)f) oc «;(0)(xi - X2) (19) 

The value of the anomaly defines the limiting contribution of the viscous term in the 
steady state: 

\im^vu{x)u" {x) = —ao{x) 

An interesting feature of this relation is that it defines the expectation value of the Oq 
- anomaly self-consistently from the steady state equation. This feature is preserved for 
the higher anomalies of the u"' - densities. They are necessarily non-zero, because after 
point splitting procedure we get terms~ -§^{u{x + f) ~ '"(^ ~ 2))*^ steady state 

equation will determine their value. 

With a certain amount of vulgarization one can say that the reason for the w" anomalies 
is that shock waves absorb not only energy, but these higher densities as well. 

Before computing the general anomaly let us discuss carefully all the limiting pro- 
cedures involved. As we see correlation functions depend on the parameters k(0) = e, 
L, (which defines the correlation length of the forces) and the viscosity u. We made an 
assumption that as we let the viscosity go to zero correlation functions have a finite limit 
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and hence depend only on e and L . This hmit is what is meant by the inertial range. The 
standard Kolmogorov assumption (which we don't make) is that the Gahlean - invariant 
correlation functions, such as {{u{xi) — u{x2))''^) have a finite limit as L ^ oo. As we 
will see in our case this statement doesn't hold. Instead we have to make a different 
assumption, consistent with our equations. We will call it the G- (Galilean) assumption. 

To formulate it let us notice that Galilean invariance in our system is spontaneously 
broken. This is evident from the fact that pumping forces create a certain average velocity 
Vrms = \/ (w^). At the same time unbroken G-symmetry would require that the probability 
distribution be invariant under u=^ u-\- const. 

It is easy to estimate the value of Vrms- In the limit of zero viscosity the only possible 
formula is : 

Vrms ^ i^^y^jL^ e^L3 (20) 
In the A representation this breakdown means that: 

EA.--^ (21) 

'-'rms 

This absence of G-symmetry makes the anomaly computation difficult. Fortunately it is 
consistent with eq.(3) to assume that if we formally tend L ^ oo and keep Xj finite, then 
G-symmetry is restored. We conjecture that in this limit: 

Z(AiXi; . . . Ajva;jv) oc 5(^Afc) (22) 

This is the G-assumption and a short check shows that it is consistent with eq.(3). To 
state this assumption in a slightly more physical way, we can say that the probability 
distributions of velocities, W{uiXi . . . u^xn), which are Fourier transforms of the Z func- 
tions, have the following structure: 

W{uiXi . . . unXn) = w{ui - Uj] Xk)Wi C^^\ (23) 



provided that : 

I Ui - Uj |< ^Uk (24) 

This last condition is very important. It is easy to see that without it the separation of 
the center of mass velocity which occured in (23) would contradict eq.(5). 

The G-symmetry greatly simplifies computations of the anomaly. However one more 
self-consistent assumption is needed . This is an assumption of the existenceof an oper- 
ator product expansion or the fusion rules. To formulate it we introduce the following 
notations: 

Z{..) = {ex^{xi) ...exj,{xN)) 

e\{x) — ex^ \u{x) (25) 
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The fusion rules is the statement concerning the behavior of correlation functions when 
some subset of points are put close together. We conjecture that in our case the rules 
have the form: 

y y d 

+ ~ ^"^^ y)(^\i+X2{x) + A2, y)^eAi+A2 + 0{y^) (26) 

We will call this statement the F-conjecture. Here A and B are some functions to be 
determined and the meaning of (26) is that they control the fusion of the functions Z^- 
into functions Z^-i as we fuse a couple of points together. To find the result one has 
to substitute (26) into (25). Of course one must check that this conjecture is consistent 
with the eq. (5), which is also is supposed to determine functions A and B. To make this 
equation effective we have to evaluate the following anomaly operator: 

ax{x) — lirn i/(A'u"(a:) exp Am(x)) (27) 

which appears on the right hand side. As we explained above, a\{x) is generally non-zero, 
because smallness of u is compensated by the blow-up of u"{x + y) exp Xu{x) as y — > 0. 
In fact we can write: 



- A ^ hm ^ ^e,(x + y)e,{x) (28) 

and exploit the F-conjecture to evaluate the RHS of (28). Thus, ax{x) should be expressed 
in terms of derivatives of the functions A and B. If the result is finite it must have the 
form: 



ax{x) = a{X)ex{x) + p{X)^ex{x) (29) 

which is the only possible G-invariant expression, involving the ultraviolet finite operators 
e{x) and u'{x). In order to have a finite limit in (28) one has to set the cut- off values 
of y and ( to be dependent on u and A. This is not surprising, since j defines the local 
Reynolds number and the ultraviolet cut-off must depend on it. It is also worth stressing 
again that this form of the anomaly is correct only in the Galilean- invariant limit. For 
generic A^ we would obtain a superposition of exponents with different A , a rather difficult 
situation to treat. 

The master equation (5) now takes the form (for the steady state): 

HZ ^ - - E - ^o)\\Z = Y: a{X,)Z (30) 

(3{X) = p{X) + J 
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It is now a simple matter to check our F-conjecture. If we introduce the variables: 

A = Ai + A2 
11 = Xi- X2 

and keep y much smaller than all other distances, we find the following structure of the 
operator H : 

Hn — Hm-i + F 

^ ^ - ^^^^'^ - ^^^'^^i + ^^^^^ " l^^^^'^ + ^^^^^^^ ^^^^ 

Here the operator Hj^_i is obtained from Hj^ by replacing the points xi 2 with the point 
X and Ai 2 with A = Ai + A2. From this we derive equations for the functions A and B 
that appear in (26): 



(Vi - V2)^ = -iPiXi) + /3(A2)) - /?(Ai + A2) 
oy 2 

dA 

(Vi - Vs)-^ = a(Ai) + a(A2) - a{X, + A2), (32) 

Here: 

These equations have solutions for any functions a and /?. The next step is to sub- 
stitute these solutions back in eqs. (28) and (29) and try to find constraints on a and 
(3. Surprisingly the arising constraints are very weak due to the possibility to adjust 
the cut-off's, and thus the functions a and /? in (29) remain almost arbitrary. In princi- 
ple they must be determined from the conditions that all probability distributions have 
admissible behaviour at 00, much in the same way in which eigenvalues are usually de- 
termined. Since we don't have any general methods for treating this problem we will 
simplify the matter even more by introducing a scaling conjecture ( S-conjecture) which 
again turns out to be self- consistent. To formulate it let us notice that if ,t ^ L we can 
expand k,{x) = k{0){1 — J2) and in (30) the constant part , k{0) drops out due to the 
G-invariance: 

^ K{xi - xj)\Xj = k{0) (^ A,)' + 0{x^) 

J2K = (33) 

Therefore it is natural to look for a scaling solution with A ~ ^. The scaling condition 
determines the possible form of the functions a and (3. In order to conform to scaling , 
they must be as following: 
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/3(A) = I 
a{X) — a 

We will see now that scaling is self-consistent, although one can also try more general 
solutions, say with logarithmic terms. 

Let us see how the unknown numbers a and b are determined from the eigenvalue 
problem. For example consider a 2-point function. The master equation in this case takes 
the form: 

For reasons to be clarified later we are interested in the case a — 0. Our S-conjecture 
amounts in the anzats: 

Ziii,y)^^Hy) (35) 

Here we temporarily use the units in which k(0) = 1 and L = 1 The function ^(x) satisfies 
an ordinary differential equation: 

x^"(x) + (1 - 2b)^'(x) + x'^^x) = (36) 

The general solution of this equation has the form: 

^x)^x^R^{^x^ (37) 

where Fx is one of the Bessel functions. The right function and the value of b are deter- 
mined from the condition that the probability distribution: 



w(u,y)= J^Z{fi,y)e-^^ (38) 

Jc-ioo ZTTl 



must be positive and vanish as w — > ±oo. From the convergence of (38) it follows that we 
must chose 

2 

(^{x) <xx''K2t{--x^^) (39) 
Positivity of Z for x > and its finitness ai x — forces us to take 6 = | for which case: 

2 3 

$(x) = exp -X2 (40) 

As a result we obtain the following result for the probability to have a velocity difference 
u at the distance y : 

w{u,y)= [ ^ewi^ifJ'y)^ - f^u] (41) 
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The positivity of w is guaranteed by the fact that Z satisfies certain convexity conditions. 
Indeed, from the relation: 

/oo 
^i;(M)e^" > 
-oo 

we derive , by using Cauchy inequahty: 

Z(/.i)^(/i2) > Z\t^^) (42) 



which is clearly satisfied. It is mathematically curious that the u representation the eq.(34) 

4 



is easily reduced to the Schroedinger equation with the potential V{u) = ^ — 2bu. This 



is achieved by the change: 

Z{fi,y) = ifi)''Fifi,y) 

, which removes the ^ term in the eq. (34) and by the Fourier transform to the function 
F{u,y) = y'^^~^4>{^)- The power of y here is needed for consistency with (35). The value 
of b corresponds to the zero energy eigenvalue in this potential and 4> is proportional to 
the ground state wave function Q. 

The probability distribution defined by (41)has the following asymptotic behaviour: 

{1 (• u \3 
e-^(y) if H ^ 
3 _^ .rl (43) 
y2u 2 II y ~^ ~00 

This qualitatively fits the observations [5 ] . It must be stressed however that some 
caution is needed when comparing the G-invariant part of the probability distributions 
with the experiment. As we have already said, the factorization (24) breaks down at large 
velocities. That means in particular that in general w{u,y) has the following structure: 

uL u y 
w{u,y) = x{ ; ; y) 

The scaling limit, discussed above is reached only when two conditions are satisfied: 

y <^ L;U <^ Vrms 

When computing the moments of the probability distributions, which represent correlation 
functions, it may happen that even when the first condition is enforced, the second one 
will be violated. In this case the result is not universal, since the behaviour of probabilities 
at ~ Vrms depends on the correlations of the stirring forces at x ~ L. 

We come to the conclusion that the breakdown of Galilean invariance leads to a rather 
peculiar structure of the correlation functions. They contain in general both universal and 
non-universal parts. The former comes from the distribution (43) and its generalization for 



^A. Migdal informed me that he was able to find this wave function directly from the Schroedinger 
equation 



10 



an arbitrary number of points. The latter results from the region of the large velocities. 
These nonuniversal correlations are just the "condensates" introduced in ref.[l]. The 
formula (43) shows that due to the "power tail" all the expectation values of starting 
with = I are formally divergent. That simply means that they are dominated by the 
non-universal region and thus change if we change k{x) for x ^ L. At the same time, the 
moments with n < | are universal. The "power tail" in the formula (43) must be related 
to the probability of having a kink, introduced in ref.[5]. However the precise connection 
is not completely clear, since the non-universal part may be relevant in the comparison. 

We come to the conclusion that at least in the present setting the violation of the 
naive scaling for higher moments - the phenomenon usualy called "intermittency"- is 
due to the breakdown of Galilean symmetry and non-universality of the large velocity 
fluctuations. In the past "intermittency" essentially meant that the theory sometimes 
works and sometimes doesn't. Here we have it under control. This observation explains 
an apparent discrepancy between the scaling in the eq. (41) and Kolmogorov's relation 
(19). These two come from the different regions of the phase space. This is evident from 
the fact that the value of k{Q) which enters into Kolmogorov's relation (19) simply drops 
off in the G-invariant limit, as seen from (33). For the different type of the stirring forces 
considered in [5] the two regions seem to overlap. That forms the basis for the beautiful 
physical picture advocated in [5]. It is also consistent with the approximate solution of 
the Burgers problem found by the replica method in [6]. 

In the above solution we took a = 0. Our understanding of other possible solutions 
is still incomplete, although it seems that for the considered type of stirring forces an 
attempt to take a 7^ leads to some unphysical results, like having {u{xi) — u{x2)) 7^ 0. 
However for a different type of stirring forces which lead to different scaling laws we almost 
certainly have to include the a term. This question is currently under investigation [ 8]. 

Finally, let us present the generalization of the master equation for arbitrary dimen- 
sionality. In the inviscid limit it isn't difficult. Consider the following quantity: 

ex = p(x,t)e^^>'*) 

where p and v are the density and the velocity, satisfying the Euler equations: 

P + da{pVa) = 

piVa + {Vi3di3)Va) = fa 

It is straightforward to verify by the same methods which led to eqs. (5) and (9) that the 
correlation function of the 's : 

F = F{Xi,xi; . . . Ajvfjv) 

satisfies the following equation: 
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which generahzes eq.(8) for an arbitrary dimension. As in one dimension the origin of 
this equation hes in the special conservation laws analogous to (13). In general we have 
the following set of conserved tensors: 



which satisfy a continuity equation. 

The next step should be an analyses of anomalies along the same lines as above. This 
task is not completed yet. Another immediate problem is to include the pressure as a 
small perturbation. This is possible to do by using the relation: 



This relation allows us to express the perturbations of pressure and density in terms of 
the function F. However this analyses is also a problem for the future. 

I am deeply grateful to V. Borue for useful discussions, A. Migdal and V.Yakhot for 
sharing with me their insights, results and enthusiasm concerning turbulence, and to D. 
Gross for his important critical remarks on physics and style of this article. 
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